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ABSTRACT 

An identical pair of thermosphere probes measured the 
Ng concentration and temperature, the ion composition and 
concentrations and the electron temperature up to 290 kilometers 
about 30 minutes and 5 minutes before totality during the March 
7, 1970 eclipse at Wallops Island, The rockets travelled sim- 
ilar trajectories thus permitting the purelv temporal changes 
between flights to be resolved. The neutral temperature and 
N 2 concentration changed little but the electron tempera- 
ture decreased by as much as 20% in the lower F-region, The 
ion concentration decreased by about 30% in the F-region and 
about 50% in the E-region, with little change in relative ion 
composition. The electron cooling rates decreased by a factor 


# 


*'%■ 

M 


r 


t 


iii 


of 6 in the lower F-region, approximately in proportion to 
the change in visible solar disc. A smaller than expected 
decrease in the cooling rate below 150 kilometers between 
the two flights indicates a hardening of the solar spectrum 
and suggests a significant heat contribution from the solar 
corona near totality. The ion composition measurements 
were consistent with solutions of the ion continuity equations. 
A proper fit required a factor of three enhancement of 
the flux below 200 A, an amount also consistent with the elec- 
tron heat balance analysis. Reactions involving the minor 

4 * + 

ions N and were found to be important for the ion 

+ + 

chemistry of the major ions O and NO , especially at 
the time of eclipse. The negligible response of the neutral 
atmosphere to the eclipse is reasonable considering the 
long time constants for the conductive and convective trans- 
port processes and the local nature of the disturbance. 
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INTRODUCTION 


The occurrence of a solar eclipse at Wallops Island on March 7, 1970 
presented a unique opportunity to investigate the response of the thermosphere 
and ionosphere to rapidly changing solar illumination at a low zenith angle. The 
response times of this region depend inherently upon the ion reaction rates, the 
electron, ion and neutral cooling rates, and perhaps on dynamic effects 
associated with the rapid changes occurring during the eclipse. To investigate 
these changes, two identical Thermosphere Probes, NASA 18. 104 and 18, 105, 
were launched southeastward from Wallops Island into the path of the approach- 
ing cone of totality. At the apogee altitude of 290 kilometers, the percent 
obscuration had reached 42% and 84%, respectively. The two rockets were 
launched along essentially identical paths to permit the purely temporal, eclipse- 
induced changes to be separated from any spatial variations that might other- 
wise have been encountered. 

In this paper, the resulting measurements of ion composition, total 
ion concentration, Nj, molecular nitrogen concentration n(N£) and tempera- 
ture, Tn 2> the electron temperature, T e , are presented and these data 
are analyzed in terms of the electron heat balance and ion photochemistry. 

THE INSTRUMENTS 

The Thermosphere Probe (TP) is an ejectable, vacuum-tight package 
that can carry several instruments to altitudes of about 290 kilometers when 
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launched by a Nike-Tomahawk vehicle. The TP has been described by Spencer, 
et al. (1965, 1969). Briefly, the Thermosphere Probe is a long cylindrical pack- 
age that is ejected from a clam-shell nose cone at about 70 kilometers, ascends 
to apcgee and returns to the Earth. After ejection, the TP tumbles end over end 
in a nearly vertical plane with a period of typically three to four seconds, thus 
exposing its instrument sensors to essentially all angles of attack during each 
tumble. 

The instruments employed in each of these two flights were an 
omegatron mass spectrometer, a Bennett ion mass spectrometer, and a cylin- 
drical electrostatic probe. These instruments have been described elsewhere 
(Spencer, et al., 1965, 1969), (Pharo, et al. , 1971), (Taylor, etal., 1963, 

1965). The omegatrons were tuned to N 2 and measured the height profiles of 
N 2 over the altitude range between 140 to 290 kilometers. The neutral gas tem- 
peratures were derived from the scale height of the N 2 profiles. The ion 
spectrometer (Pharo and Scott, 1971) measures the ion concentrations in the 
mass range from 12 to 36 AMU, sweeping at a rate of about 5 times per second. 

The cylindrical probe experiment employed two redundant rhodium 
plated collectors and sweep voltage rates of 22 and 80 volts per sec. A sweep 
repetition rate of 8 per second was employed. The electron temperatures, T , * 

V 

were derived from both probes in the usual manner of fitting an exponential to 

% 

the retarding region (Brace, et al 1970), The ion concentrations, Nj, were derived 
from the ion currents collected when the probe was perpendicular to the velocity 
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vector, following the method described by Taylor, et al, (1963). Details of 
the data processing for these flights will not be discussed in this paper. 

THE MEASUREMENTS 

The pre-totality launch times of 1800 UT (1300 EST) for NASA 18. 104 
and 1827 UT (1327 EST) for NASA 18. 105 were selected to permit resolvable 
changes to occur in the various measured parameters without encountering the 
non-linear rates of change that might be expected near totality (Stubbe, 1970). 
This was expected to simplify the theoretical analysis of the measured profiles. 
The area of unobscured solar disc, as viewed from the rocket trajectories, 
changed by a factor greater than four between the two flights thus assuring 
significant changes in solar energy input. As noted earlier, the obscuration at 
apogee of the two flights was 42% and 84%, respectively. Figure 1 shows the 
path of totality as it passed Wallops Island at 18:38 UT and the locations of 
apogee and the 100 kilometer ascent and descent crossings of the two flights. 

The first rocket reached apogee at 18:04 UT and the second at 18:31 UT. 

Figure 2 shows the percent obscuration along each leg of the trajectory. The de- 
scent data were selected for detailed analysis because of the generally smaller 
angles of attack for the N 2 and ion measurements and the larger percent obscur- 
ation. Figure 3 shows the n(N 2 ) profiles measured at these times and the corre- 
sponding T n profiles derived from the N 2 scale heights. The T e profiles are 
also shown. There was essentially no resolvable change in either n(N 2 ) or Tj^ 2 
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Figure 2. The Percent Obscuration Encountered Along the Trajectories 
of the Flights. Data from the Descent Leg of Both Flights 
Have Been Employed in This Paper. 
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within the 5% relative accuracy of these measurements. T c decreased at all 
altitudes with the greatest decrease, 20%, occurring at 200 kilometers. 

The ion concentration profiles from the two flights, shown in Figure 
4, were derived by normalizing the constituent ion currents to the total ion 
concentration Nj derived from the probe experiment. The latter agreed well with 
the N e profiles derived from simultaneous ionosonde records taken at Wallops Sta- 
tion (Jackson and McQuillan, 1970) except in the region where the Nj values were 
10-15% lower. The 5-10% agreement in the E and F^ regions is considered ac- 
ceptable considering the horizontal distance between Wallops and the descent 
legs of the trajectories. Tables I and II are tabulations of the measurements 
taken on descent of the two flights. 

ACCURACY OF THE MEASUREMENTS 

The factors that affect the absolute and relative accuracy of these 
measurements have been treated elsewhere, but it is appropriate to reiterate 
those factors that are especially important for the eclipse study. Since the 
goal of this investigation was to resolve the changes in each parameter during 
the eclipse we have attempted to maintain high relative accuracy by making 
the two sets of instruments as identical as possible. For this reason the 
omegatrons were calibrated simultaneously for on the ground pressure 
calibration system normally used for this purpose. Their relative accuracy 
should be better than 5%, including the possible errors in knowledge of the angle 


of attack as measured by earth and sun sensors on both flights. The T N 

2 

accuracy is also believed to be better than 5%. 

The accuracy of the probe measurements of T g were estimated to be 
5% on the basis of the quality of the exponential fits to the retarding regions 
and the good agreement between the measurements from the two identical probes 
on each flight. The Nj values were estimated to be accurate within 10% based 
on their agreement with the ionosonde N @ profiles In the ascent E- region and 
F-region essentially above Wallops Island. 

The concentrations of the constituent ions are expected to be ac- 
curate within 20%. Before they were normalized to the Nj profiles the indi- 
vidual ion currents were corrected for mass discrimination within the spec- 
trometer using the same factors as used in previous applications (Pharo, 
et al., 1971). 

The Neutral Atmosphere Model Adopted 

Since the physical interpretation of these data relies in part on a 
knowledge of the neutral composition of the thermosphere at the times of these 
flights, we have adopted the Jacchia Model that most closely matches the 
n(Ng) and T^ profiles. Figure 5 and Table m represent the model used in 
the subsequent analyses of these data. To aid in matching the model, we 
extends! the N 2 and T^ 2 profiles downward by interpolation with the 
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ION CONCENTRATION (ions/cm 3 ) 

Figure 4. The Ion Composition Profiles. 
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Figure 5. The Jacchia Model Atmosphere Best Matching the 
n(N2) and Profiles Measured on These Flights. 
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essentially simultaneous Pitot-probe measurements of Horvath and Theon (1971) 
who measured the total atmospheric density and temperature up to 120 kilometers. 

Magnetic Disturbance 

It should be noted that the eclipse occurred during a period of large 
magnetic disturbance. Table IV lists the 3-hour magnetic index, ap, for a three 
day period centered on the eclipse day. ap was 07 at the time of these flights 
(see box in Table IV). This could be an important factor in the interpretation 
of the eclipse data, as there is increasing evidence that storms affect the tem- 
perature and composition of the neutral atmosphere (Jacchia, et al, 1967) at 
ionosphere heights thus altering the ion chemistry and thermal balance. We 
have not attempted to consider the storm effects explicitly, however, 

DISCUSSION 

The response of the Thermosphere to a solar eclipse is of considerable 
interest as the solar obscuration introduces abrupt and, to son extent, calcul- 
able changes in the energy input to this region. Observations of the thermospheric 
response to the eclipse in turn provide the opportunity to test our theoretical 
understanding of the processes controlling its energy and chemical balance. 

These have been examined in the past using data obtained under normal daytime 
conditions (Brace, et al, 1969) (Pharo, et al, 1971). In this paper we employ 
the same type of analysis to eclipse time measurements. In particular we 
have examined primarily the changes in concentration and temperature observed 
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between flights to evaluate the changing energy input and photochemical 
balance during the onset of the eclipse. 

The Energy Balance 

The magnitude of the electron temperature in the lower ionosphere is 
primarily determined by the local equilibrium between the electron heating and 
cooling rates (Hanson, 1962) (Dalgarno, etal,, 1963), The electrons are 
created as photoelectrons that, in their thermalization process, heat the ambient 
electron population. The latter are cooled by elastic and inelastic collisions 
with ions and neutrals. Although the heating rate is difficult to calculate lacking 
simultaneous solar euv information, its equivalent, the cooling rate, can be 
calculated readily from the measurements of these flights. In this section, we 
consider the electron heat balance in an attempt to evaluate the effects of the 
changing solar obscuration during this eclipse. 

Following Brace, etal., (1969), the T 0 , T N2 » Nj, and n(N 2 ) profiles, 
and the selected Jacchia model atmosphere (Figure 5), were employed to calcu- 
late the electron cooling rates during these flights. We included the effects of 
heat conduction (Banks, 1966) and the ion heating induced by thermospheric 
winds (btubbe and Chandra, 1970), with a range of wind velocities assumed 
independent of height. The temporal variation of the internal energy was 
neglected, as a simple calculation shows that it contributed less than 1% to the 
energy budget for the small changes that occurred during the eclipse. 


Figure 6 shows the resulting cooling rates (104) and Q e (105) from 
the two flights* where the family of curves represents the range of wind veloci- 
ties assumed. The effects of winds become increasingly important above 220 km 
where the winds may cause Ti to exceed Tg. The effects of heat conduction 
represent an additional uncertainty in the cooling rates at high altitudes. The 
heat conduction flux depends strongly upon the height gradient of T e , a factor 
that c annot be derived accurately at the higher altitudes where the horizontal 
motion of the rocket and temporal changes may be comparable to the altitudinal 
variation o* T e . However, below 220 kilometers, the wind and heat conduction 
effects become negligible compared to local cooling, and the validity of the cool- 
ing rate profiles depends only upon the accuracy of the measurements and the 
model. 

Figure 6 also compares the ratio of the cooling rates from the two 
flights (solid) with the ratio of the unobscured solar disc areas (dashed). If the 
heat sources arose entirely from radiation from the visible disc, these ratios 
would be comparable, and this appears to be approximately true in the range 
from 150 to 220 kilometers. However, at lower altitudes, the cooling rates did 
not decrease proportionally to the decrease in the visible solar disc, suggesting 
that some of the energy for this region originates outside the visible disc. This 
is not surprising because solar X-rays originating in the corona are a known 
source of E-region ionization and heating. Similar evidence for a '’residual 
flux" was reported by Smith, et al., (1965) who found little change in T e at lower 
altitudes during the South American eclipse of 1967. 
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Figure 6. The Cooling Rotes Calculated from NASA 18.104 and 18.105 
Profiles and Their Ratio Q e (104)/GL (105). Wind Heating of 
the Ions Is Included for Speeds of 0, TOO, 200 Meters /Sec. 

The Ratio of Unobscured Solar Areas A, / A« Is Shown for 
Comparison. 
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The Ion Chemistry 

In our analysis of the ion chemistry, we adopted the ultraviolet spectrum 
of Hitteregger, et al., (1965) and scaled it upward linearly in accordance with 
the higher level of solar activity at the time of the eclipse (F^ 0 ^ - 171). We 
then included the solar obscuration effect at each point along the trajectory and 
exposed the model atmosphere to these euv fluxes and, following Pharo, et al. 
(1971), calculated the height profiles of each ion constituent, This 
approach was modified to consider possible transient effects by including in the 
continuity equations the changes observed between flights in the various ion 
constituents. This correction was particularly important for the second flight 
where it accounted for about 30% of the 0 + concentration above 200 kilometers. 
For all other ions, the effect was less than 4%. From this it is perhaps valid 
to generalize that a steady state solution of the continuity equation is appropriate 
even during an eclipse for all ions except 0 + . 

The neutral atmosphere model used for the ion chemistry analysis was 
identical to that used in the cooling rate calculations, and the T e measurements 
were employed to calculate the temperature dependent recombination rates. 

These initial calculated ion profiles (not shown) were found to agree 
with the measured ion profiles within a factor of two. This might be considered 
reasonable agreement in view of the uncertainty in the magnitude of the euv flux, 
the hardening of the euv spectrum implied by the cooling rate profiles, and the 
uncertainty in the absorption and ionization rates. The most pronounced 
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deficiency in the initial theoretical result was the low 0 + concentration below 
200 kilometers, a deficiency that in turn affected the other ion concentrations. 

To attain a better fit to the measurements, two modifications were 
made. First, the absorption and ionization rates were increased by a factor 
of two from the values given by Hinteregger et al. (1965) and used by Pharo 
et al, (1971) on non-eclipse measurements. These changes are within their 
quoted uncertainties. Secondly, we decreased the rate coefficient for the 
reactions 

n + + N 2 — ■ NO + + O 

from 

6. 5 x 10" 13 to 3 x 10" 13 cm 3 sec” 1 

and 

o + + o 2 — o 2 + + O 

from 

5.9 x 10"’ 12 to 3.5 x 10” 12 cm 3 sec” 1 . 

Both of these rates are lower than laboratory derived rates (Dunkin et al., 1968) 
by abr :t a factor of two. 

To introduce the amount of hard radiation that was inferred from the 

smaller than expected change in cooling rates, the euv spectrum was scaled up- 

« 

ward by a factor of three below 200 A, as shown schematically in Figure 7. 

Since the visible solar disc was only 7% in the E-region of the second flight, 
this factor of three scaling suggests that 14% of the radiation at these wavelengths 
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Figure 7. Schematic Model Used for Additional Hard Radiation (Matched) 
Needed to Account for the Excess Electron Heating Rates and 

Deficiency of O on the Second Flight. A Factor of Three 
Enhancement Was Assumed at Wavelengths Below 200 A. 
Detailed Spectrum Is Not Shown. 
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originated outside the obscured area. From this one would predict that the 
x-rays are only about HCft obscured by the moon during totality. 

The effect upon the O* population of intensifying the hard radiation by 
this scheme is shown in Figure 8, The relative enhancement of n(O f ) is shown 
for three altitudes assuming the upper wavelengths of 100, 200, 300, 400, and 


e 

500 A, respectively. 


The wavelengths above 300 A are relatively less effective 


at low altitudes. Most of the desired enhancement results from wavelengths 
between 100 and 200 A. 

The final theoretical ion profiles are shown (dashed) in Figures 9 and 
10, and are compared with the measured ion profiles (solid). The agreement 
is generally better than 15% and thus appears to fall within the combined 
experimental uncertainties. 


Uniqueness of the Analysis 

It should be recognized that the analyses that lead to the theoretical 
ion profiles oi Figures 9 and 10 is not entirely unique with respect to the 
absolute values of ion concentration. Perhaps the greatest uncertainty lies 
in our ignorance of the solar euv spectrum at the time of the eclipse. Had we 
used a spectrum different from that presented by Hinteregger et al (1965), 
somewhat different adjustments of the ionization, absorption, and reaction 
rates would have been necessary to fit the observations. 

However our conclusion that a hardening of the spectrum occurred 
during the eclipse is on firmer ground as it is based solely on the relative 
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RELATIVE 0* ENHANCEMENT 



UPPER LIMIT OF ENHANCEMENT (A) 


Figure 8. The Relative Enhancement of 0 + at 135/ 170 and 240 Kilometer 
Caused by the X-ray Enhancements of the Type Shown in Figure 

, 7 Assuming Upper Limits of 100, 200, 300, 400 and 500 A. Th<, 

J Major 0 + Enhancement at Low Altitudes Occurs for Wavelengths 

f Between 100 and 200 A. 
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Figure 9. Comparison of 18.104 bn Measurements (Solid) with 
Theoretical Profiles (Dashed) Calculated for the 
Corresponding Solar Obscuration, Including the 
Enhanced X-ray Component. 
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changes in the ion concentrations and electron temperatures between the two 
flights. The changes in these parameters require enhancements in the 100- 
200 A range. The form we chose to introduce this enhancement (Figure 7), is 
not unique, however. 
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The Effects of th e Minor Ion Reactions 

The calculations of Stubbe (1970) for this eclipse disagree significantly 

with regard to the behavior of ion composition during the eclipse. Although 

some of the discrepancy could be accounted for by suitably modifying the solar 

euv input, one reason for the failure in Stubbe’ s model to predict the NO + and 

Or? correctly is because it does not consider the role of the minor constituents 

N + and N*> + in the following reactions: 

N + + 0 2 — NO + + O 
N + + o 2 — *0 2 + + N 

V + °2“~°2 + + N 2 

N 2 + + O— NO + + N 

Both of the N + reactions are important as the N + concentrations are 
significant in the lower ionosphere and their rates are of the order of 5x 10“^® 
cm 3 sec" 1 (Fehsenfeld, et al., 1965) (Ferguson, 1967). To illustrate the impor- 
tance of N + , we have included in Figure 10 the solutions in which these reactions 
have been neglected (dotted). Without N + , the Og + concentrations are a factor 
of 3 too low at 160 km and NO + is low by almost a factor of 2. The N 2 + reactions 
have less influence on the O ^ and NO + concentrations but are not negligible. 
They account for about 20% of the 0 2 + and NO + . Pharo, et. al., (1971) have 
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Figure 10. Comparison of 18.105 Ion Measurements with Theoretical Profiles, 

with X-rays Included. The Dotted Lines Represent the O 2 + and NO + 
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included these reactions in their analysis of the normal daytime thermosphere 
above Wallops Island and obtained similarly good fits to the composition 
measurements. 

The Neutral Atmosphere 

It is perhaps appropriate to comment upon the apparent lack of response 
of the neutral atmosphere to the eclipse. As noted earlier, no resolvable 
changes in n(N 2 ) or occurred between these flights. However, owing to 
the magnetic disturbances occurring at the time of the eclipse (a p - 67) and 
throughout at least the previous 24 hours, one should perhaps not assume 
that the response to the eclipse was entirely normal. Any energy transport 
from higher latitudes would tend to mask the eclipse effect and further contribute 
to a smaller response of the atmosphere. 

For two other reasons it remains understandable that the eclipse 
caused little change in the neutral atmosphere. The first is related to the long 

response time of the neutral atmosphere and the second is related to the highly 
localized nature of the eclipse. 

The total internal energy in a 1cm 2 column above 120 km is about 
Q«5x 10 4 ergs cm” 2 . 

This energy can be conducted downward at the 120 km level with a flux of 2.7 x 
10“ 1 ergs cm" 2 sec" 1 with a resulting characteristic time of about 1.85 x 10 5 
sec or about 51 hours, a time much greater than the 27 minutes between the 
flights. Furthermore, vertical and horizontal heat convection (Volland and 
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Mayr, 1970) can supply energy to and from the adjacent atmosphere at rates 
of 5 x 10" 1 erg cm" 2 sec" 1 and 4 x 10“ 1 erg cm" 2 sec” 1 , respectively, leading 
to characteristic times of the order of 1 x 10® sec. These also are much longer 
than the eclipse time, thus it is not surprising to find little change during the 
eclipse. Since the solar obscuration affects only a very small region at any 
given instant, horizontal transport tends to fill up this "energy hole" from the 
adjacent and fully illuminated regions without significantly affecting the global 
energy budget. The effect of this process is to further damp the eclipse 
response of the thermosphere. This aspect has been considered by Volland and 
Mayr (1971) who show that small scale energy disturbances are ineffective in ex- 
citing density variations. The excitation efficiency varies as l/n 2 , where n is 
the wavenumber corresponding approximately to the ratio of the global dimen- 
sion to the size of the disturbance. 

Stubbe (1970) predicted that the exospheric temperature would decrease 
during the eclipse, at the end of the eclipse falling about 120° lower than its 
normal day value. However, Stubbe’s model overestimates the variations to 
be expected during the eclipse. It is one dimensional and therefore does not 
consider the damping effects of horizontal mass and energy transport, factors 
that are expected to be significant. 
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CONCLUSIONS 


Measurements of Ion, electron and neutral particles at two times dur- 
ing the eclipse have permitted the electron heat balance, the ion chemistry and 
the neutral atmospheric response to be evaluated. Both the electron heat balance 
and ion chemistry suggest a hardening of the solar spectrum as the eclipse 
progresses, with X-rays attenuated about a factor of three less than the euv 
flux. The euv flux appears to have been attenuated in proportion to the visual 

obscuration. The role of the minor ions N and N g is found to be important 

+ + 

in the production of NO and O , These minor ir: is are usually neglected, but 

are especially important during an eclipse when photoionization sources for 
+ + 

NO and O are less significant. The lack of variation of neutral temperature 

A 

and concentration appears consistent with the long thermal and transport time 


constants of the thermosphere. 
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Table III 




Jacchia Model. 

T* “ 1120°K 
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Table IV 

a„, 3- Hour Magnetic Index 


B333S 

0-3 

3-6 

6-9 

9-12 

12-15 

15-18 

18-21 

21-24 

March 6 

27 

15 

15 

22 

7 

18 

27 

67 

7 

56 

32 

15 

27 

32 

39 

■K31I 

67 

8 

32 

48 

67 

48 

154 

207 

400 
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